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Abstract 

In supersymmetric quantum mechanics, exact-solvability of one-dimensional quantum systems 
can be classified only with an additional assumption of integrability, the so-called shape invariance 
condition. In this paper we show that in the prepotential approach we proposed previously, shape 
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I. INTRODUCTION 



It is generally known that exactly solvable systems are very rare in any branch of physics. 
Thus any new method to construct exactly solvable models would be of interest to the 
community concerned. It is therefore very interesting to realize that most exactly solvable 
one-dimensional quantum systems can be obtained in the framework of supersymmetric 
quantum mechanics (SUSYQM) However, in SUSYQM, exact-solvability can be 

classified only with an additional assumption of integrability, so called shape invariance (SI) 
condition Hence in SUSYQM the SI condition must be taken as a sufficient condition for 
integrability at the outset. What is more, the transformation of the original coordinate, say 
X, to a new one z = z{x) needed in solving the SI condition is not naturally determined within 
the framework of SUSYQM in most cases, but have to be taken as given from the known 
solutions of the respective models. It would be more satisfactory if the exact-solvability of 
a quantal system, including the required change of coordinates, could be determined with 
the sirnplest, and the most natural requirements. 

In y, la, Isl a unified approach to both the exactly and quasi-exactly solvable systems 



i§, presented. This is a simp 
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e constructive approach, based on the so-called prepotential 
isl . which gives the potential as well as the eigenf unctions 
and eigenvalues simultaneously. The novel feature of the approach is that both exact and 
quasi-exact solvabilities can be solely classified by two integers, the degrees of two poly- 
nomials which determine the change of variable and the zero-th order prepotential. Hence 
this approach treats both quasi-exact and exact solvabilities on the same footing, and it 
provides a simple way to determine the required change of coordinates z{x). All the well- 
known exactly solvable models given in [l|, I2I, most quasi-exactly solvable models discussed 
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and some new quasi-exactly solvable ones (also for non-Hermitian 
Hamiltonians) , can be generated by appropriately choosing the two polynomials. 

Since all the well-known one-dimensional exactly solvable models obtained in SUSYQM, 
by taking SI condition as a sufficient condition, can also be derived without the SI condition 
in the prepotential approach, one wonders what role the SI condition plays in the latter 
approach. In this paper we would like to show that the SI condition is only a necessary con- 
dition in the prepotential approach to exactly solvable systems. Therefore, unlike SUSUQM, 



shape invariance needs not be assumed in the prepotential approach. 
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This paper is organized as follows. In Sect. II we give a brief review of the prepotential 
approach to exactly solvable models with both sinusoidal and non-sinusoidal coordinates. 
The idea of SI as a sufficient condition of integrability in SUSYQM is sketched in Sect. III. 

Sect. IV and V then demonstrate that in the prepotential approach for models with sinusoidal 
and non-sinusoidal coordinates, SI is automatically satisfied and needs not be imposed. 
Sect. VI concludes the paper. 



II. PREPOTENTIAL APPROACH 

The main ideas of the prepotential approach can be summarized as follows (we adopt the 
unit system in which h and the mass m of the particle are such that H — 2m = 1). Consider 
a wave function (pN^x) {N: non- negative integer) which is defined as 

0^(x) = e-^°(^W(^), (1) 

with 

f 1, A^ = 0; 

\Uk=ii^-Zk), N>0. 
Here z = z{x) is some real function of the basic variable x, Wo{x) is a regular function of 
z{x), and z^s are the roots of Pn{z)- The variable x is defined on the full line, half-line, 
or finite interval, as dictated by the choice of z{x). The function Pn{z) is a polynomial in 
an {N + l)-dimensional Hilbert space with the basis (1, z^, . . . , z^). Wo{x) defines the 
ground state wave function. 

The wave function (p^ can be recast as 

(f)N = exp {-WNix, {zk})) , (3) 

with given by 

N 

Wm{x, {zk}) = Wo{x) - J] In \z{x) - Zk\. (4) 

fe=i 

Operating on by the operator —(P/dx"^ results in a Schrodinger equation ifjv0iv = 0, 
where 

Hn--^,+V^. (5) 

= w'^ - w;,. (6) 



Here prime represents differentiation with respect to x. It is seen that the potential Vn is 
defined by Wn, and we shall call Wn the A^th order prepotential. From Eq. (jl]), one finds 
that Vn has the form Vat = Vq + AVn: 



^ ^ k=l " hd 



{z - Zk){z - Zi)' 



Thus the form of Vat, and consequently its solvability, are determined by the choice of 
Wo{x) and z''^ (or equivalently by z" = {dz'^ /dz)/2). Let W^^z' = P^{z) and z'^ = 
be polynomials of degree m and n in z, respectively. In Q], it was shown that if the degree 
of Wqz' is no higher than one (m < 1), and the degree of z''^ no higher than two {n < 2), 
then in Vn{x) the parameter N and the roots z^s, which satisfy the so-called Bethe ansatz 
equations (BAE) to make the potential analytic, will only appear in an additive constant and 
not in any term involving powers of z. Such system is then exactly solvable. If the degree 
of one of the two polynomials exceeds the corresponding upper limit, the resulted system is 
quasi-exactly solvable. The transformed coordinates z{x) such that the degree of z'"^ is no 
higher than two are called sinusoidal coordinates. There are six types of one-dimensional 
exactly solvable models which are based on such coordinates, namely, the shifted-oscillator, 
three-dimensional oscillator, Morse, Scarf type I and II, and generalized Poschl- Teller models 
as list in 

In jsi, the prepotential approach to exactly solvable systems was extended to systems 
based on non-sinusoidal transformed variable z{x) which is a solution of z' = X — z"^. With 
this, the remaining four types of exactly solvable systems listed in pj], namely, the Coulomb, 
Eckart, and Rosen-Morse type I and II models, are also covered by the prepotential approach. 

A. Sinusoidal coordinates 

For exactly solvable models with sinusoidal coordinates we take m = 1 and n = 2, i.e., 
Pi{z) = az + b, and Q2{z) = az^ + [3z + '~^, where a, 6, a, /5 and 7 are real constants. With 
these choices we obtain |4| 
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■ (8) 



Demanding the residues at z^s vanish gives the set of Bethe ansatz equations 

■ ^-§-5;*M.o, * = l,2,...,iV. (9) 

4 Zk — Zi 



With this set of roots Zk, the last term in Eq. ([8]) vanishes, and we obtain a potential 
Vn{x) = Vo{x) — En without simple poles. Here Vo{x) = Wq — Wq does not involve and 
ZkS, and can be taken as the exactly solvable potential of the system with eigen-energies 



En = 2aN — aN"^ . In fact, Vo(x) is exactly the supersymmetric form presented in [ij] for the 
shifted-oscillator, three-dimensional oscillator, Morse, Scarf type I and II, and generalized 
^oschl- Teller models (for easy comparison, we note that a and a here equal and aA in 

B. Non-sinusoidal coordinates 

As mentioned before, the Coulomb, Eckart, and Rosen-Morse type I and II models involve 



a change of coordinates of the_ 
extension of the methods in 



brm z' = X — z"^ which is non-sinusoidal. But with a slight 
^, all these four models can be treated in a unified way in 



the prepotential approach j^. The extension is simply to allow the coefficients in Wq be 
dependent on N. It turns out that Wq takes the form 

W;,iN) = -(A + Na)z + ^fj^, (10) 

where A and B are real parameters. Then the potential Vn becomes Vn{x) = V{x) — En, 
where 

V{x) = A{A-l)z\x) -2Bz{x), (11) 

and 

En = -jxf]:^^ - ^ ^2Ar + 1) + N'] . (12) 

Now V{x) is independent of A^, and can be taken to be the potential of an exactly solvable 
system, with eigenvalues En (N = 0, 1,2, . . .). The corresponding wave functions 0Ar are 
given by ([1]): 

0^ ^e(^+'^)^'''"^(")-^>^(a;), AT = 0, 1, . . . (13) 



The BAE satisfied by the roots z^'s are 



E 



X 



Zk - zi 



{A + N-l)zk + 
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0, k 



1,2,. ..,iV. 



(14) 



Finally,we mention here that V{x) in (fTT!) can be obtained, up to an additive constant, 
from Wq{N) with anyvalue of A^. Particularly, the form adopted in supersymmetric quantum 
mechanics (e.g., in [l|) is obtained from the zero-th order prepotential Wq{N = 0) with 



N = 



a. 



III. SHAPE INVARIANCE IN SUPERSYMMETRIC QUANTUM MECHANICS 

From the discussions in the last section, we see that in the prepotential approach, exactly 
solvable models are determined by the zero-th order prepotential Wo{x) in the sinusoidal 
cases, or Wq = Wo{N = 0) with = in the four non-sinusoidal cases. The potential 
Vo is completely determined by Wq: Vq = — Wq, and consequently, the Hamiltonian 
Hq = —dP/dx^ + Vo is factorizable as Hq = A with the first-order operators 

This fact is indeed the base of SUSYQM. In SUSYQM one considers the relation 

between the spectrum of and that of its so-called super-partner Hamiltonian Hi con- 
structed according to H\ = AA^ = —d'^/dx'^ + Vi, where Vi = Wq + Wq. In forming Vi, 
it is equivalent to using a prepotential —Wq. The ground state of Hi is therefore exp(Wo), 
and it follows that the ground states of Hq and Hi cannot be both normalizable. 

Let us suppose the ground state of Hq, i.e. exp(— l^o); is normalizable, and denote 
the normalized eigenf unctions of the Hamiltonians Hq i by ipn'^^ with eigenvalues En'^\ 
respectively. Here the subscript n = 0,1,2,... denotes the number of nodes of the wave 
function. It is easily proved that Vq and Vi have the same energy spectrum except for the 
ground state of Vq with E^^ = 0, which has no corresponding level for Vi jll, Q]. More 
explicitly, we have the following supersymmetric relations: 



^« = (i^it)-^wt, ^v^r=o, (16) 
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Hence A annihilates and converts an eigenfunction of an excited state of Hq into an 
eigenfunction of Hi with the same energy, but with one less number of nodes, while ^4+ 
does the reverse. Consequently, if the spectrum of one system is exactly known, so is the 
spectrum of the other. 

This is, however, all that supersymmetry says about the two partner potentials. If any 
one of the spectra is unknown, then supersymmetry is useless in solving them. It is therefore 
gratifying that most of the well-known one- dimensional exactly solvable models process a 
property called shape invariance. With hindsight, one can then impose shape invariance 
as an additional requirement along with supersymmetry to classify exactly solvable systems 
having such property. This has been done and most exactly solvable systems are then unified 
within the framework of SUSYQM [l, Q. 

Shape invariance means that the two super-partner potentials Vq and Vi are related by 
the relation 



Vl(x;Ao) = Vo(x;Ai)+/?(Ao), (17) 

where Aq is a set of parameters of the original Vq, Ai = /(Aq) is a function of Aq, and -R(Ao) 
is a constant which depends only Aq. This implies 

W;,\x, Ao) + <(x, Ao) = W^\x, Ai) - <(a;, Ai) + R{Xo). (18) 

Eq. ( fTTIl implies that Vi has the same shape as that of Vq, but is defined by parameters 
Ai instead of Aq. From (fTSll one deduces that the ground state wave function of Vi is 
~ exp{—Wo{x, Ai) with energy i?o(Ao). Then from (fT6|) we know the energy of the first 
excited state of Vq to be -R(Ao), and the wave function tpf^ ~ A'^ip^^^ By repeated use of the 
shape invariance condition, one can construct the partner V2 of Vi, V3 of V2, etc. The ground 
state wave function of (n = 0, 1, . . .) is V'o"'' ~ exp(— l^o(a^j A„), where A„ = /"(Aq), with 
energy Ylt=o ^i^k)- Then again from (JT6l) we know that the wave function of the n*^ state 
of Hq is ijjn^ ~ {A'^)"-ipl^\ with energy 

n-l 

i?W = ^i?(Afc), n = 0,l,... (19) 

k=0 

So with shape invariance one obtains the complete spectrum of Hq. 

It is now obvious that SI is a sufficient condition of integrability in SUSYQM. To classify 
shape-invariant exactly solvable models in SUSYQM, one must solve the SI condition (fTSj) 



to get all the functional forms of Wo{x), Ai = /(Aq), and i?(Ao). This general problem 
is very difficult and, to the best of our knowledge, is still unsolved. Further constraints 
on the possible class of shape invariant potentials are required. Particularly, in order to 
obtain the well-known exactly solvable models one must assume that (again with hindsight) 
the parameters of the two partner potentials are related by simply a translational shift, 
i.e. Ai = /(Ao) = Ao + m differ from Aq only by a set of constants m. Even with this 
simplification, the required change of coordinates z = z{x) needed in solving the SI condition 
cannot be determined naturally in the approach of SUSYQM, but has to be taken as given 
from the known solutions of the respective models. 

On the other hand, in the prepotential approach SI needs not be imposed, and Wq and 
z{x) are determined by simply picking two polynomials with the appropriate degrees. In 
this sense it appears to us that the prepotential approach is conceptually much simpler. 
Nevertheless, putting the differences of the two approaches aside, one could not help but 
wonder what role SI plays in the prepotential approach. Below we would like to demonstrate 
that for the exactly solvable models obtained in the prepotential approach, SI is automat- 
ically satisfied. We shall discuss the cases with sinusoidal and non-sinusoidal coordinates 
separately. 

IV. SHAPE INVARIANCE IN PREPOTENTIAL APPROACH: SINUSOIDAL 
COORDINATES 

Our strategy is to show that, with z{x) and Wo{x) given in Sect. 11(A) and (B) that 
produce the ten well-known exactly solvable models, the SI condition f|T8l) is always satisfied, 
i.e. one can always find the set of new parameters Ai in terms of the old ones Aq. In the 
process, we demonstrate that the change in the parameters of the shape-invariant potentials 
are translational. 

In this section, we first consider the cases involving sinusoidal coordinates. For exactly 
solvable systems, we must take Wqz' = Pi{z). LabeUing the corresponding parameters of 
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the two shape- invariant potentials by = 0, 1, we have 

z'^ = Q2{z) = az^ + (3z + r, (20) 
Pi'\z) = akZ + h, , k = 0,l (21) 

W^{Xk) = ^=^, \u = {a,,h,). (22) 

yQ2{z) 



Note that z[x) is the same for the shape-invariant potentials. Then the SI condition (ITSll 
leads to 

{pr - pr) + (pf' + Pi'') - Y-^ {pp + Pi'') = Ri^o)Q. (23) 

Equating the coefficients of the powers of z, one arrives at the following equations relating 
the parameters 

al — a1 = Ra, 

2 {aobo - ai6i) + ^ (ao + ai) - a {bo + 6i) = Rp, (24) 
6^-62 + 7 (ao + ai) - ^ (&o + &i) = ^7- 

For simplicity we write R for R{\q). We mention here that the signs of a and b are fixed 
by the normalization of the wave functions. This means they are the same for the two 
shape-invariant partner potentials. 

We would like to solve (12^ for Ai = (ai,6i) and R in terms of Aq = (ao,&o)- To 
facilitate solution, we find it convenient to first determine all inequivalent types of sinusoidal 
coordinates. 



A. Inequivalent sinusoidal coordinates 

Depending on the presence of the parameters a, (3 and 7, there are three inequivalent 
cases of sinusoidal coordinates: (i) 2;'^ = 7 7^ 0, (ii) z'"^ = Pz + •y {(3 7^ 0), and (iii) 
z''^ = az'^ + /3z -|- 7 (a 7^ 0). By an appropriate shifting and/or scaling, these cases can be 
recast into three canonical forms. 

The form given for case (i) is already the canonical form of this case. We shall take 7 > 
as 7 < leads to physically uninteresting change of variable. This case gives rise to the 
shifted oscillator. 



9 



By shifting z to z = z + 'y/f3 in case (ii), we get the canonical form z''^ = f3z. For physical 
systems we require /3 > 0. This case corresponds to the three-dimensional oscillator. 

Case (iii) can be recast as 5'^ = azP' + 7, where z = z + (3/2a and 7 = A/4a with the 
discriminant A = 4q;7 — For the case A = (the exponential case) and a > 0, the 
system thus generated is related to the Morse potential. For A 7^ 0, we have two situations. 
If a > (the hyperbolic case), the canonical form is z'"^ = a(i^ ± 1), where z = >/4a^J\A\z, 
and the plus (minus) sign corresponds to A > (A < 0). The plus sign gives rise to the 
Scarf II model, while the minus sign corresponds to the generalized Poschl- Teller model. 
For a < (the trigonometric case), the canonical form is z''^ = |a|(±l — i^), where again 
z = A/4a^/|A|z, and the plus (minus) sign corresponding to A < (A > 0). With the plus 
sign we get the Scarf I model, while the minus sign does not lead to any viable system as 
the transformation is imaginary. 

From the above discussions, we see that we need only to discuss the three inequivalent 
canonical cases, namely, (i) z'^ = 7 7^ 0, (ii) z'"^ = Pz {P > 0), and (iii) z''^ = a{z'^ + 6) 
{5 = 0, ±1 for a > 0, and 5 = -1 if a < 0). 

B. Case (i): z'2 = 7 > 

For this case, it is easy to check that ao (o-i) must not vanish, or it will lead to vanishing 
potential. Furthermore, we must have ao > and ai > in order that the wave functions 
be normalizable. The SI conditions ( 12^ become 



Equations (125|) and (126|1 require Oi = Oq, &i = bo, or ai = — Oq, &i = — &o- In the latter 
solution the signs of ai and bi are different from those of ao and b^, and hence the wave 
functions of one of the two systems cannot be normalizable if those of the other system can. 
In fact, for this case we have R = from ( 1271) . This means the ground states of the two 
systems have the same energy. But the flip of both signs of a and b of Wq means that the 
ground states of the two systems have the forms exp{—Wo) and exp(+Wo)- They cannot be 
both normalizable. This is exactly the result in SUSYQM. 



^0-^1+7 (oo + oi) = -R7. 



(ao + ai)(ao — ai) = 0, 



ao^o Ciibi = 0, 



(25) 
(26) 
(27) 
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So we are left with the choice oi = oq, bi = Bq. From (127|) we have R = 2aQ. Thus R is 
a constant, and from f|T9l) it imphes oscillator-hke spectrum, i.e. = naQ. This gives the 
shifted oscillator. 

The above discussion shows that in this case SI is a necessary condition. The parameters 
of the two partner systems are related by (oi, bi) = (oq, 60); and the shift parameter is 
R = 2ao. 



C. Case (ii): z'^ = pz {p > 0) 

Normalizability of wave functions in this case require that a > and b < 0. Now the SI 
conditions (HM are 

(ao + ai)(ao - ai) = 0, (28) 
2 (ao&o - fli&i) + ^ (ao + ai) = (29) 
(60 + &1) (^o-&i-0 =0. (30) 

Possible solutions of these equations are ao ± ai = 0, 60 + ^1 = or 60 ~ ^1 ~ /3/2 = 0. To 
keep the signs of a and b unchanged, we can only take (ai, 61) = (oq, bo — P/2) as the viable 
solution. Then from fl30|) we get R = 2ao, which again gives an oscillator- like spectrum. 
This is just the case of the three-dimensional oscillator. 



D. Case (iii): z'^ = a(z'^ + S) 

Next we consider the case with z'^ = a{z'^ + 5) [S = 0, ±1 for a > 0, and 5 = —1 if 
a < 0). As mentioned before, this case covers the Morse, generalized Poschl- Teller, and the 
Scarf I and II potentials. The SI conditions (|2^ are 

al — ai = Ra, (31) 
2(ao6o-ai6i)-a(6o + &i) = 0, (32) 
bl ~ bl + aS (flo + «i) = RdS. (33) 

To solve ai , bi and R in terms of Qq and bo, we eliminate Ra in ([3SD using fl^ to get 
{bo + bl) {bo - bl) + 5 {ao + ai) {ai - ao + a) = 0. (34) 
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From (IMIl we can have four possible sets of solutions: 

ao + ai = 0, bo + bi = 0; (35) 

aQ + ai = 0, bo — bi = 0; (36) 

oo — cti = a, bo + bi = 0; (37) 

ao — cti = «, 6o ^ ^1 = 0. (38) 

The first three sets of solutions involve change of signs of a and/or b, and so are not viable 
as discussed before. Thus for this case we must take (ai, bi) = (ao — a, 6o) which also 
satisfies (15^ . Eq. (PTi) then gives 

i?(Ao) = = 2ao - a. (39) 

From (fT9!) the energies are 

O-o 



En 



a 

Oq — (ao — na)'^ 



n = 0, 1, . . . (40) 



a 

This is exactly the results in SUSYQM 

To conclude this section, we have shown that SI is automatically satisfied in the prepo- 
tential approach for the sinusoidal cases. 

V. SHAPE INVARIANCE IN PREPOTENTIAL APPROACH: NON- 
SINUSOIDAL COORDINATES 

In this case, Wq = —Az + BjA and z' = a{X — z"^). Here Ao = {A,B). As in the last 
section, we show that one can always find a set of new parameter Ai = (A', B') in terms of 
Ao that solves the SI condition f|T8|) . In fact, from f|T8|) one finds 

A{A + a) = A' {A' -a), (41) 

B = B\ (42) 
B^ B'^ 

— - aXA = ^ + OiXA' + R. (43) 

Solutions of (HTj) are A' = —A and A' = A + a. The first solution has the sign of A 
changed, and will lead to non- normalized wave functions. Hence the viable solution is 
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Ai = {A',B') = {A + a, B). Once again, the change in the parameters A and B of the 
shape-invariant potentials are translational. Finally, from fH3|) we find 



B' 



a\{2A + a). 



(44) 



_A^ {A + a)\ 

This agrees with the results in SUSYQM 

Thus we have shown that in the prepotential approach for models based on non-sinusoidal 
coordinates, SI is also a necessary consequence of the forms of Wq and z'. 



VI. SUMMARY 



A unified approach to both the exactly and quasi-exactly solvable systems has been pro- 
posed previously based on the so-called prepotential in 4,5,^. In this approach solvability 
of a quantal system can be solely classified by two integers, the degrees of two polyno- 
mials which determine the change of variable and the zero-th order prepotential. All the 
well-known exactly solvable models obtained in SUSYQM can be easily constructed by ap- 
propriately choosing the two polynomials. 

But all these exactly solvable models are obtained in SUSYQM only by taking the SI 
condition as a sufficient condition. The requirement to get exactly solvable models in the 
prepotential approach appears to be much simpler, and definitely without the need of SI 
condition. In this paper we have shown that the SI condition is in fact only a necessary 
condition in the prepotential approach to exactly solvable systems, and hence needs not 
be assumed. In the process, we have demonstrated that the change in the parameters of 
the well-known shape-invariant potentials are indeed translational, a result which was also 
assumed in SUSYQM. 
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